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Abstract 

Shape-phase transitions in mixed parity boson systems are studied using mean 
J> \ field theory. It is shown that without parity projection, shape transitions in mixed 

parity systems are very similar to those in positive parity systems, both requiring a 
finite interaction strength for the onset of deformation. The shape-phase diagram in 
0"N mixed parity systems, however, changes significantly after parity projection, with 

the critical point for the onset of octupole deformation moving to zero strength. 
Shape transitions in the Lipkin model are shown to exhibit the same features, 
indicating that this is a direct consequence of parity projection in finite, mixed 
parity systems independent of the nature of interacting particles. 
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The low-lying negative parity states in even-even collective nuclei are described 
in terms of the octupole degree of freedom (see [1,2] for recent reviews). Ex- 
perimental study of the octupole states is hindered by the much stronger 
quadrupole collectivity which limits their accessibility. Thanks to the advances 
in detector technology, this situation has changed dramatically in recent years. 
Exploiting the capabilities of the 7-ray arrays Eurogam and Gammasphere, 
the range of negative parity spectra in actinide nuclei Rn, Ra and Th has been 
greatly extended [3,4]. These new data, in turn, provide fresh challenges for 
collective models, the most interesting one being the description of the con- 
jectured transition from the octupole-vibrational to the octupole-deformed 
shapes in actinides. In this letter, we present mean field studies of this shape- 
phase transition in two analytically solvable models, the interacting boson 
model (IBM) [5] and the Lipkin model [6], that will shed further light into 
this question. 
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We first briefly review the spherical — > quadrupole-deformed shape transition 
in the sci-IBM [7-9], which provides a useful reference point for the octupole 
case. A convenient form of the Hamiltonian is 

H sd = e d h d - k 2 Q ■ Q, (1) 



where the one-body term tries to keep the system spherical while the quadrupole 
interaction drives it towards deformation. Here h d = dj^d^ is the rf-boson 
number operator and Q = [s^d + d^s] + Xd[d^ d]^ is the quadrupole operator. 
For simplicity, we use the axially symmetric intrinsic state in the mean field 
analysis (the full state leads to the same result) 

\N,P 2 ) = 



2\N 



s f + M |0>. 



(2) 



Evaluating the expectation value of the Hamiltonian (1) in the intrinsic state 
(2) yields the energy surface 



E(fo) = N 2 k 2 



rj 2 P; 



i + Pi 



' 2(3 2 + XdPl 



(3) 



where Xd = — y 2/7xd, and r\ 2 = e d /NK 2 provides a convenient parameter to 
study the shape transition. Here and in the following, the 1/7V corrections are 
neglected in the spirit of the large N assumption inherent in the mean field 
theory. Variation of E(f3 2 ) with respect to (3 2 gives 



fa mix + Pi) - 2(2 + xM(i + Xdfo - Pi 



o. 



(4) 



Analysis of Eqs. (3, 4) shows that the critical point for spherical — > deformed 
shape transition is given by i] 2 = 4 + x% F° r Vz > ^ ne system has a 
spherical shape with j3 2 = 0, and for r\ 2 < rj 2 , it becomes deformed with (5 2 
determined from Eq. (4). At the critical point, the spherical and deformed 
solutions coexist with energies E(f3 2 = 0) = E((3 2 = Xd/2) = 0. We note that 
angular momentum projection leads to the same energy expression as (3) to 
leading order in iV [10], hence restoration of the rotational invariance leads to 
a qualitatively similar shape-phase diagram. 

The proper framework for discussion of shape-phase transitions in mixed par- 
ity systems is the spdf-IBM which has been shown to give an adequate de- 
scription of the positive and negative parity spectra in collective nuclei [11-14]. 
However, the mean field analysis of the spdf-IBM is quite complicated which 
tends to obscure otherwise a very simple conclusion. To illustrate the main 



2 



ideas involved in a simple case, we first discuss the spherical — > octupole- 
deformed shape transition in the s/-IBM using the Hamiltonian 

H sf = £ f hf - k 3 • O, (5) 



where n,f is the /-boson number and O = [s^ f '+ p s]® is the octupole operator. 
This Hamiltonian is similar to (1) when Xd = 0, and like the 7-instability 
in that case, it is unstable in the direction of deformation. That is, each of 
the octupole mean fields [3 3fl is populated equally, and the energy surface is 
a function of j3 2 = With this provision, the intrinsic state and the 

corresponding energy surface (without PP) are given by 



\N,p 3 ) = [N\(l+Pl) N 
E((3 3 ) = N 2 k 3 



-1/2 



(V + zW^lo), 



2ft 



l + PI \i + p 2 3/ 



(6) 



where 773 = €f/NK 3 . Variation of E{(3 3 ) w.r.t. f3 3 gives 

/3 3 ^3(1 + - 4(1 - = 0, 



(7) 



which has the solutions 



1) /3s = 0, E = 0, 

2) Pl = ^ £ = ~^3-4) 2 , , 3 <4. (8) 

4 + ^3 Id 

It is clear from these solutions that for r\ > 4, the shape is spherical, and 
the critical point for the onset of octupole deformation occurs at r] 3 = 4. 
Thus, just like in the quadrupole finite interaction strength, given by 

K 3 — e //4-V, is required for the shape transition to take place. Note that 
one can obtain this result without explicit solutions by simply expanding the 
energy surface around fi 3 , and looking at the coefficient of f3 3 . From Eq. (6), 
E((3 3 ) ~ N 2 k 3 [(t] 3 — A)f3 3 + ...], and it is obvious from the sign of r] 3 — 4 that 
the (3 3 = solution is stable for r/ 3 > 4 but becomes unstable for r/ 3 < 4, 
leading to a deformed solution. This method of determining the critical point 
is useful in cases where the extremum conditions cannot be solved explicitly. 

The intrinsic state in Eq. (6) has mixed parity, and it is important to see 
how PP influences the shape transition characteristics discussed above. The 
projection operator onto good parity is given by 

P n = (l + nV)/2, (9) 
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where tt = ±1 is the parity quantum number and V is the parity operator. 
Under V, boson operators transform as Vb\ m V = (—l) l b lm , hence V acting on 
the intrinsic state (6) gives 

V\N,f3 3 ) = \N,-(3 3 ). (10) 

From Eqs. (9) and (10), the projected states can be written as 

P v \N,fo) = (\N,p 3 )+n\N,-p 3 ))/2. (11) 



Binomial expanding the even and odd spin parts in (11), it is easy to see that 
for 7i = +1, it has an even number of / bosons, and for ix = — 1, it has an odd 
number of / bosons. Normalization of the projected state (11) is given by 

■Mr = (N,[3 3 \P w \N,f3 3 ) = (l + nr N )/2, (12) 



where 

r = (l-(3 2 3 )/(l + (3 2 3 ), (13) 



gives a convenient measure of parity mixture in the condensate; for r = 1 
{(5 3 = 0), it consists of only s bosons, while for r = {(5 3 = 1), it is an equal 
admixture of s and / bosons. 

The energy surface of the Hamiltonian (5) with PP is given by 



E n (p 3 ) = ^ r (N,P 3 \HP n \N,p 3 ), 



Variation of Eq. (14) leads to a complicated equation that cannot be solved 
analytically. Nevertheless, by plotting the energy surface for various values 
of 773, one can gain useful insights into the nature of the shape transition. In 
Fig. 1 we show plots of Eqs. (6) and (14) for r] 3 = 4, 40 and 400. The first value 
corresponds to the critical case without PP. The next two values represent a 
10- and 100-fold reduction in the octupole strength compared to the first case 
for a fixed €f/N. In the top figure, the energy surface without PP (dashed 
line) exhibits the typical flat-bottomed curve expected at the critical point. 
After PP the energy surface for i] 3 = 4 is seen to become deformed (solid line). 
What is surprising, however, is that despite the orders of magnitude reduction 
in the octupole strength in the middle and bottom figures, the f3 3 = solution 
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remains unstable and the minimum of E + always occurs at j3 3 > 0, that is, 
the energy surface after PP remains deformed. As k 3 is progressively reduced, 
the well depth and deformation decrease proportionally but the energy surface 
itself (solid line) stays more or less shape invariant. To find the critical point 
analytically, we expand E + (14) around /3 3 = 



E + ~ N 2 k 3 -2$j + (773 - 2)Pl + ... 



(15) 



It is clear from Eq. (15) that for any k 3 > 0, there is a deformed minimum, and 
hence the critical point is given by k 3 = 0. That is, the onset of the octupole 
deformation is immediate in the presence of any octupole strength, and the 
spherical phase found in the study of shape transition without PP disappears 
after PP. This situation is depicted in Fig. 2, where /3 3 corresponding to the 
absolute minimum of the energy surface with and without PP are plotted 
against I/773. The sharp spherical — > octupole-deformed transition exemplified 
by the dashed curve (which is also representative of the quadrupole case) is 
replaced by a smoothly varying curve after PP, and the critical point moves 
from I/773 = 0.25 to (or from k 3 = e f /AN to 0). 

To understand this somewhat surprising result, it is worthwhile to look at the 
effect of PP on individual terms by comparing Eqs. (6) and (15). The leading 
order (/3f) contribution to the one-body energy is seen to vanish after PP, 
whereas the octupole term survives, loosing only half of its strength. A more 
direct way to see the effect of PP is to expand the condensate 

+(1/2)N(N - l)(sY- 2 0s ■ P) 2 + • • • , (16) 

and look at the overlaps of various terms. Denoting the states in the expansion 
by their / boson number rif, the leading contribution to the one-body term 
comes from the (1-1) matrix element (m.e.) (i.e. which goes like 

However, after PP, odd nf terms are projected out, so the leading term has to 
come from the (2-2) m.e., which goes like f3 3 . In contrast, the leading (equal) 
contributions to the octupole term come from the (1-1) and (0-2 + 2-0) m.e., 
and PP blocks only the first one, leaving the second one intact. Thus the net 
effect of PP is to lift the obstruction of the one-body operator against the 
formation of a deformed system. 

The above conclusion on the effect of PP on shape transitions in mixed parity 
systems can be generalized to arbitrary boson spins, and in particular to the 
spdf -IBM. However, as stressed before this case is rather complicated and 
is deferred to a more detailed study. Due to weaker dipole collectivity, the p 
boson is weakly coupled to the sdf system, and should only have a perturbative 
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effect on the sdf -IBM results. Thus, for our purposes, it is sufficient to study 
the physically most relevant case, namely, the onset of octupole deformation 
in a quadrupole deformed system within the framework of the sdf '-IBM. The 
Hamiltonian is given by a combination of Eqs. (1) and (5) 

H sdf = e d h d + Efhf - k 2 Q ■ Q - K3O.O, (17) 



where the quadrupole and octupole operators are generalized to 
Q = [ a tJ +(i t a ](2) +Xd [Sdf) + X /[/ t /] (2) , 

= [ S t/+ /t s ](3) + X3[d 1f + /tjj(3). (18) 

Numerical studies of this Hamiltonian with a general intrinsic state show that 
the ground state remains axially symmetric for realistic choices of parameter 
[15]. Using the intrinsic state 

\N,(3 2 ,(3 3 ) = -± ; 1/2 |0), (19) 

the energy surface of the Hamiltonian (17) after PP is given by 



"(2ft + Xd (3l + XfPl? + *r N -\2{3 2 + Xd fi - xs 0lf 



K>2 



-^(2[3 3 + 2 X3 {3 2 p 3 ) 2 j, (20) 

where (5% = 1 + $ + r = (1 + /3| - /5|)/^, Xf = ^A/21 Xf and xs = 

— ^4/15x3- The energy surface without PP follows from Eq. (20) by setting 
7r = 0. In the case of no PP, an analysis similar to Eqs. (7, 8) can be carried 
out leading to an octupole vibrational phase (^3 = and f3 2 given by Eq. (4)), 
and a deformed phase (^3 > 0) which is too complicated to present here. 
Nevertheless the critical point can be obtained from the expansion of the 
energy surface around f3 3 = as 

_ e/(l + PI) + 2NkM2 + Xd f3 2 ){{3 2 -x f ) (0U 

K3 ~ m + X3&) 2 • (21) 

The quadrupole deformation is seen to have modified the critical point for 
the onset of octupole deformation but still a finite k 3 is required. A similar 
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expansion of Eq. (20) for it = +1 shows that the leading term is proportional 
to — k 3 /?| as in the spherical case (15). Thus after PP, the critical point again 
occurs at k 3 = 0. 

As a final example, we consider shape transitions in the Lipkin model [6] , which 
has been shown to mimic, at the mean field level, the octupole deformation 
effects obtained from the more realistic Hartree-Fock-Bogolyubov calculations 
[16]. In the Lipkin model, iV = 2 j + 1 nucleons occupy two levels that have 
the same angular momentum j but opposite parity n. The Hamiltonian with 
one- and two-body terms is given by 

H = eK - ^k(K + K + + K_K_), (22) 

where X's are the SU(2) quasi-spin operators defined as 
1 „ / * 



— 7} Z_, \ a m+ a m+ — a m- a m-) , 
m 

K+=J2 a ln+a m -, K_ = {K+)\ 

m 

with the quantum numbers k = N/2, ji ■ 

In the "spherical" (or harmonic) limit (k 
are given by 

i0o>=n«uo>, 

m 



(23) 

= -N/2, -N/2 + 1,..., N/2. 
= 0), the ground state and its energy 

E Q = -eN/2. (24) 



The excited states are obtained by repeated application of K + on the ground 
state, which generates a finite, harmonic spectrum, E n = E + ne, n = 
1,...,N. For k > 0, the two body-interaction mixes the n = +1 and — 1 
levels, and to find the ground state of the system, one needs to use an intrin- 
sic state with mixed parity 



\N,(3}= (l + (3 2 ) 



j2\N 



-1/2 



II («L + Pal + ) |0>, 



(25) 



where the mean field (3 describes the "deformation" of the system. With the 
intrinsic state (25), the energy surface of the Hamiltonian (22) without PP is 
given by 



E{P) = ~ + N\ 



vP 2 



P 



l + (3 2 \l + (3 2 



(26) 
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where 77 = e/Nn. Apart from the constant term, the energy surface (26) is 
very similar to that of the s/-IBM (6), and it is easy to see from its expansion 
that the critical point for the onset of deformation occurs at rj = 1 or k = e/N. 

Parity projection in the Lipkin model is defined just like in the boson case 
Eqs. (9-11) and yields the following energy surface 



EM 



l + P 2 



1 + ixr 



N-2 



1 + ixr 



N 



(27) 



where r is defined as in Eq. (13). Expanding E + {j3) around (3 — gives 

E + ~ —eN/2 — N 2 k(3 2 + Thus the critical point again moves to k = 

after PP. 

In conclusion, using analytically solvable models of mixed parity systems we 
have demonstrated that a sharp transition to octupole deformation, similar to 
the onset of quadrupole deformation, occurs only if one ignores PP. After PP, 
the character of the shape transition changes completely, and the vibrational 
phase disappears. While one can still use the terminology and formalism of 
octupole vibrations for small fa, this is only an approximation and not a well 
defined phase as in quadrupole vibrations. A similar effect of the restoration of 
discrete symmetries in finite systems is well known in the BCS theory, where 
the critical pairing strength vanishes after the number projection. The effect of 
PP on shape transitions, however, seems to have been ignored in the literature 
so far. Implications of this result for the octupole states in actinides will be 
discussed in a longer article. 

Parts of this work were carried out while visiting the Japan Atomic Energy 
Research Institute and the Institute for Nuclear Theory at the University 
of Washington. Hospitality and support of these Institutes are acknowledged 
with thanks. 
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Figure Captions 

Fig. 1. Energy surface with ir = +1 (solid line) and without PP (dashed line) for 
various values of 773 ranging from 4 to 400. The former corresponds to the critical 
value without PP. The energy surface with n = — 1 is also shown in the top figure 
(dotted line) but it is too high in energy to show in the other cases. 

Fig. 2. Shape-phase diagram for spherical — > octupole deformation with PP (solid 
line) and without PP (dashed line). Deformations, obtained from the minima of the 
energy surfaces Eqs. (6) and (14), are plotted against 1/7/3 = Nns/ef. 
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